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QUANTIFYING THE ANISOTROPY IN THE INFRARED EMISSION OF POWERFUL AGN 
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ABSTRACT 

We use restframe near- and mid-IR data of an isotropically selected sample of quasars and radio 
galaxies at 1.0 < z < 1.4, which have been published previously, to study the wavelength-dependent 
anisotropy of the IR emission. For that we build average SEDs of the quasar subsample {— type 1 
AGN) and radio galaxies (= type 2 AGN) from ~ 1 — 17 /zm and plot the ratio of both average samples. 
From 2 to 8 /im restframe wavelength the ratio gradually decreases from 20 to 2 with values around 
3 in the 10 /im silicate feature. Longward of 12 /im the ratio decreases further and shows some high 
degree of isotropy at 15 /im (ratio ~1.4). The results are consistent with upper limits derived from the 
X-ray/mid- IR correlation of local Seyfert galaxies. We find that the anisotropy in our high-luminosity 
radio-loud sample is smaller than in radio-quiet lower-luminosity AGN which may be interpreted in 
the framework of a receding torus model with luminosity-dependent obscuration properties. It is also 
shown that the relatively small degree of anisotropy is consistent with clumpy torus models. 
Subject headings: galaxies: active - galaxies: high-redshift - infrared: galaxies - galaxies: nuclei 



1. INTRODUCTION 

In the unification scheme of AGN the difference be- 
tween type 1 and type 2 AGN is explained by angle- 
dependent circumnuclear o bscuration of th e accretion 
disk and broad-line region (IAntonucci|[l993D . This ob- 
scuring dusty medium - commonly referred to as "dust 
torus" - is optically and geometrically thick and prob- 
ably extends from sub-parsec scales outward to several 
10s of parsecs, or beyond for high luminosity objects. 
The dust in the torus absorbs the incident UV/optical 
radiation an re-emits the received energy in the infrared. 

Observations have shown that type 1 AGN show sig- 
nificantly more emission in the near-IR tha n ty pe 2 A GN 
for th e same given intrinsic luminosity (e.g. lLeipski et al.l 
120051 ) . This is consistent with the picture where the face- 
on view onto the torus in type 1 AGN exposes the inner- 
most hot dust to the observer. On the other hand in type 
2 AGN the torus is seen edge-on so that internal obscu- 
ration blocks the line-of-sight to the hot dust. Owing to 
this effect, it is expected that for a given AGN luminosity 
the infrared emission of type 1 AGN is generally stronger 
than from type 2s. 

In the light of attempts at forming isotropic AGN sam- 
ples based on IR fiuxes it seems important to know ex- 
actly how strong of a bias towards type Is over type 2s 
may occur when invoking fiux limits. Moreover, probing 
the wavelength dependence of this anisotropy in the in- 
frared has some constraining power on our understanding 
of how the torus obscures the AGN. It may be possible to 
distinguish torus models where the dust is smoothly dis- 
tributed from t hose where the dust is arranged in clouds 
(e.g. Lcvenso n et al.ll2009| ): If the dust is smoothly dis- 
tributed within the torus, a large degree of anisotropy is 
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expected. If, however, the dust is arranged in clouds the 
anisotropy is expected to be smaller. 

A problem commonly encountered when studying 
AGN samples in the local universe is a significant contri- 
bution of the host galaxy to the IR. This is related to the 
typical lower luminosity Seyfert galaxies which dominate 
the nearby AGN population. One way around this prob- 
lem is the use of high-spatial resolution observations, as 
possible with the largest ground-based telescopes or in- 
terferometers, which are able to resolve out the host and 
isolate the AGN emission (for details see iHonig et al.l 
I2010D . However, it is difficult to set up representative 
samples owing to the observational limitations. Another 
possibility is the use of high luminosity objects - typi- 
cally at higher redshift - where the AGN outshines the 
host galaxy by a large factor in the optical and near- 
IR. If PAH features are absent, the AGN most likely 
dominates the mid-IR wavelength region as well (in our 
sample: host ^ 10% for wavelengths < 17 /im). 

In this paper we aim at quantifying the wavelength de- 
pendence of the anisotropy of the AGN emission in the 
infrared from ~ 1 — 17 /im. For that we use a (nearly) 
isotropically selected and complete sample of quasars and 
radio galaxies with hidden quasars at 1. < z < 1.4 as 
recently presented in iLeipski et al.l ()2010D and described 
in Sect. [2] Here we improve the analysis by using host- 
galaxy subtraction for the radio galaxies and use clumpy 
torus models for interpretation. In Sect. |3] we show the 
average SEDs of each of the subsamples which are repre- 
senting obscured (type 2) and unobscured (type 1) AGN. 
We further analyze the origin of the anisotropy by fitting 
extinction models and clumpy torus models to the obser- 
vations in Sect, m In Sect. [5] we discuss our results by 
comparing them to previous IR anisotropy estimates in 
literature. The results are summarized in Sect. [6l 

2. SAMPLE SELECTION AND DATA 

The object sam ple for this paper comprises all 3CRR 
(|Laing et al.lll983[ ) radio galaxies and quasars with 1.0 < 
z < 1.4 This lobe dominated sample presents a well 
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matched set of radio galaxies and quasars in terms of 
their intrinsic luminosity ((j^L^^i7g mhz) = (1-9 ± 0.6) x 
10^*erg/s for the quasars and {I'L^^nsuHz) — (2.0 ± 
0.4) X 10''^ erg/s for the radio galaxies; errors indicate 
standard deviation of the sample ). The data used here 
have been presented previously in lHaas et al.l ()2008t ) and 
iLeipski et al.l (|2010D and we refer the reader to these pa- 
pers for further details on the source selection, data re- 
duction, and building of the average SEDs. To sum- 
marize briefly, we obtained mid-IR photometry in all 
six filters from 3.6 /im to 24 /im and spectroscopy from 
19 /im to 38 /im utilizing all three instrurnents o nboard 
the Spitzer Space Telescope ([Werner et al.|[2003 ). After 
the data reduction in a standard manner, the individ- 
ual source SEDs were interpolated onto a common rest 
frame wavelength grid. The quasars and radio galax- 
ies were then averaged into a mean SED for each class 
of objects. The individual SEDs (including observed and 
interpolated pho tometry) as w e h as t he the average SEDs 
are presented in ILei pski et al.l (|2010l their Figs. 1 & 2). 

For this paper, additional corrections have been ap- 
plied to the radio galaxy data before the averaging pro- 
cess outlined above: At the shortest wavelengths con- 
sidered here (^1 — 3 /im, rest frame) the radio galaxy 
SEDs show contributions from the host galaxy. Since 
we want to isolate the emission coming from the active 
nucleus, we have to correct for the stellar emission in 
these cases. This correction was performed by fitting 
the observed IRAC photometry with a combination of a 
moderately old (5-10 Gyr) elliptical galaxy SED to rep- 
resent thfi^tenaremission (taken from the GRASIL web- 
page; iSilva et all I1998D and a hot black body for the 
AG N dust whose temperature we allowed to va ry (see 
e.g. iSevmour et all [2007L iDe Breuck eTaD IMol ). The 
resulting black body temperatures in the radio galaxies 
range from 600 to 970 K, with a median value of 860 K. 
The fraction of host galaxy light contributing to the flux 
measured at the observed frame wavelengths 3.6, 4.5, 5.8, 
and 8.0 /tm was found to be 0.9, 0.6, 0.3, and 0.1, respec- 
tively. Despite a negligible contribution from the host 
galaxy, we performed similar fits to the quasar sample 
as well to obtain characteristics of the hot dust emission 
and compare it to the radio galaxies. For the quasars 
we find notably hotter temperatures in the range from 
880-1250K (median 1020 K). This is consistent with the 
idea that the hottest dust is obscured in radio galaxies, 
while directly seen in quasars. 

In the radio galaxies we then subtracted the estimated 
host galaxy contribution from the observed IRAC pho- 
tometry using the scaled template SED. For the IRS and 
MIPS measurements the host galaxy contributions were 
considered negligible at restframe wavelength > 7 /im 
and no corrections have been applied. We refrain from 
further correct ions related to possib le starformation. As 
pointed out in ILeipski et al.l (j2010[ ) neither the individ- 
ual SEDs nor the averages showed any PAH features (see 
also Fig. [T]), which indicates that any contribution from 
starformation to the mid-IR is probably negligible. 

The radio galaxy 3C 368 has a galactic M-star super- 
imp osed close to the position of the radio galaxy nucleus 
(e.g. lHammer et al.lll991l:lBest et al.lll997D . Both sources 
are partly blended even at the shortest IRAC wave- 
lengths which makes the correction for the host galaxy 



emission in this source quite uncertain. Consequently, 
we removed this source from the sample considered here. 
This leaves us with 11 quasars and 8 radio galaxies from 
which the average SEDs have been calculated. 

3. INFRARED SED OF QUASARS AND RADIO 
GALAXIES AT Z - 1.2 

In Fig. [1] we show the average SED of 1.0 < z < 1.4 
quasars (red) and radio galaxies (blue) respectively. The 
error bars reflect the mean absolute deviation while the 
shaded areas show the range of the respective subsam- 
ple at each wavelength point of the interpolated data 
(see Sect. [5]). For each of the object types we calcu- 
lated a spline fit through the mean data points in order 
to guide the eye. It is obvious that the radio galaxies 
are systematically lower in infrared emission than the 
quasars. The discrepancy is largest in th e near-IR and 
flatten s out towards longer wavelengths. ILeipski et al.l 
(|2010[ ) showed similar average SEDs. Here we used addi- 
tional host galaxy subtraction, which isolates the AGN 
light much better (see Sect. [2|). This is most obvious in 
the near-IR part of the radio galaxies shortward of 5 /im. 
The SED keeps on falling toward shorter wavelengths 
consistent with the Wien tail from h ot dust emission, in- 
stead of making an upward turn (see ILeipski et al.ll2010l 
Fig. 2). 

Since both types have the same radio luminosities due 
to our selection, this difference between quasars (= type 
1 AGN) and radio galaxies (= type 2 AGN) is a generic 
property of the sample. It either reflects a difference 
in line-of-sight extinction (e.g. by cold dust in the host 
galaxy) , or traces the anisotropy in r e-emission of the 
AGN-heated dust. ILeipski et al.l (|2010[ ) tested the former 
possibility and found a surprisingly good match of the 
difference between radio galaxies and quasars in the mid- 
IR by a single extinction law. This, however, breaks 
down in the near-IR. In this paper we will test if, instead, 
a single absorber and emitter (= the dust torus) may be 
responsible for the radio galaxy /quasar anisotropy (see 
Sect. 1121). 

To quantify the wavelength dependence of the 
anisotropy we plot the quasar/radio galaxy ratio in 
Fig. [21 Also shown is the mean absolute deviation of the 
ratio calculated by propagating the standard deviations 
of each subsample. From 2 to 8 /im the emission ratio 
gradually decreases from 20 to 2. In the silicate feature 
this ratio increases again up to about 3 and flattens out 
towards 15 /im at a value of 1.44 ± 0.17. 

Our sample comprises a range of SED shapes, meaning 
that a range of ratios is observed. Most of this sample 
range comes from the radio galaxies which show much 
less uniformity than the quasars (see Fig. [JJ. We illus- 
trated the range of ratios covered by the radio galaxies in 
Fig. [3| where we plot the wavelength dependence of the 
ratio of each radio galaxy using the average quasar SED 
and normalize it to 15 /im. 

At around 10 /tm the radio galaxies show the silicate 
feature in absorption while quasars display a weak sili- 
cate emission feature. We used the spline fits described 
above (see Fig. [TJ to locate the centers of the silicate 
features following the method outlined bv iSirockv et al.l 
l|2008D . The silicate absorption feature center is found to 
be at 9.6 ±0.1 /im while the peak of the silicate emission 
feature was measured at 10.6±0.1 /tm. This difference in 
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Fig. 1. — Infrared sample average SED of quasars (red filled cir- 
cles) and radio galaxies (blue filled circles). The gray error bars 
represent the dispersion (mean absolute deviation) of each sample 
while the red- and blue-shaded areas indicate the total range cov- 
ered by the quasar and radio galaxy sub-samples, respectively. The 
red-dotted and blue-dashed lines are cubic spline fits to the quasar 
and radio galaxy average SED. 

central wavelength is similar to the ones observed in lo- 
cal galaxies. As recent high spatial resolution studies of 
Seyferts suggest, the "shift" in wavelength is not a pure 
radiative transfer effect due to the location or distribu- 
tion of the dust, but implies a chang e of dust chemistry 
wi thin the torus (iHonig at al.ll20ldf ). As demonstrated 
bv lSmith et all (pCmTT the silicate emission feature may 
be located at around 10.5 /^m if a fraction of the hot sil- 
icate dust consists of porous silicate grains. 

We note that the transition from quasars to radio 
galaxies is not as smooth as one may expect. In spite 
of some overlap in the range of SEDs of both types in 
the near-IR part in Fig. [1] quasars generally show in- 
frared emission characteristics expected for a type 1 AGN 
(blue IR color; silicate emission feature) and radio galax- 
ies have IR SEDs with type 2 emission properties (redder 
SED; silicate absorption feature). 

4. ANALYSIS 

In this section we analyze the scenarios that can lead 
to the observed anisotropy in radio galaxies and quasars. 
There are two possibilities: (1) extinction by a cold dust 
screen and (2) absorption and emission in a warm dusty 
medium. The former possibility may be associated with 
cold dust in the host galaxy while the latter one is equiv- 
alent to the dust torus in the AGN unification scheme 
(which we may call "intrinsic anisotropy" ) . 

4.1. Cold dust screen extinction 

We will first discuss the plausibility and consequences 
of a cold dust screen on the SEDs. If the IR anisotropy 
is dominated by cold host galaxy dust, then we have a 
situation where the IR emission originates from the torus 
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Fig. 2. — Wavelength dependence of the ratio of quasar to radio 
galaxy SED as a measure of anisotropy of the infrared emission. 
The black-filled circles are the ratios of the average SEDs as shown 
in Fig.[T]and the dark-gray error bars are the mean absolute devia- 
tion for each wavelength bin. The red-dashed line is a comparison 
to clumpy torus models. The dark-blue dotted line shows a stan- 
dard ISM dust screen extinction curve scaled to the type 1/type 
2 ratio at 15 fira, while the light-blue dotted lines is sc r een e xtinc- 
tion using the IR absorption curve by IChiar fc Tieleni l|2006l ) . For 
details see text. 



while the absorption is coming from a different compo- 
nent (i.e. dust in the host). This means that we require 
an additional component outside the AGN to model the 
data (e.g. as used for a sig i iifican t minority of high-z type 
2 QSOs in lPolletta et al.l ()2008[ )V However, the objects 
suffering extinction (here: radio galaxies) would be offset 
from quasars by only a sing le extinction law. T his has 
been tested and ruled out bv lLeipski et al.l (j201 Cfi for the 
same set of objects as presented here. In Fig. [2] we show 
the anisotropy ratio between quasars and radio galax- 
ies (black circles with error bars). We overplot a stan- 
dard ISM extinction curve, resembling cold screen ex- 
tinction, scaled to the observed 15 fim anisotropy (dark- 
blue dotted line; using a mixture of 53% silicates and 
47% gr a phite , based on updated dust opacity cuves by 
iDrainel 0OOl. a nd a grain size distribution according 
to iMathis et all (|1977f )). The extinction curve signifi- 
cantly overpredicts the anistropy in the silicate feature 
with better agreement in the near-IR. For reference we 
also p lot the extinction curve based on' Chiar fc Tieleni 
(120061 light-b l ue do tted line; "Pixie dust"), as used in 
iLeipski et al.l ()2010[) , which results in better agreement 
within the silica te feature but signi ficant offsets in the 
near-IR. In fact, ILei pski et al.l ()2010[ ) pointed out that a 
good correspondence of quasars and radio galaxies can 
be achieved only if the quasar SED is attenuated by at 
least two instead of one extinction components - which 
is reminiscent of radiative transfer (= absorption and 
emission) within the torus rather than a cold screen. 

Moreover, we found a mean anisotropy of about 1.4 at 
15/Lim (see Sect.[3l). Dust opacity curves typically have 
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Fig. 3. — Wavelength dependence of the ratio for each individual 
radio galaxy using the mean quasar spectrum <F(Quasar)>. All 
individual ratios have been normalized to 15 fim (and the IRS range 
> 10 fim has been binned) for illustrating the range covered by our 
sample. 

opacity ratios Tis^m/Ty ~ 0.01. In order to obtain the 
observed anisotropy, the dust screen would have to have 
Ty ~ 30. While such optical depth values are in reach 
for galactic dust lanes (e.g. extinction towards our own 
Galactic center) , it requires very edge-on views onto disk 
galaxies since the scale height of galactic disks is small. 
On statistical grounds this possibility may be viable for 
a minority of all radio galaxies, but it is unlikely that the 
whole population is dominated by host extinction. We 
note that the same line of argument can be made using 
the near-IR anisotropy leading to even higher Ty and il- 
lustrating the need for more then just one cold extinction 
screen in this scenario. 

4.2. Torus model fitting 

We use our clumpy torus models CAT3D 
(|Honig fc Kishimotol I2010( ) to test if the observed 
anisotropy ratio can be explained by the intrinsic 
anisotropy as predicted in the unification scheme. 
Generally smooth dust torus model s predict stronger 
aniso tr opy than clumpy m odels (jSchartmann et aLl 
[20081) . iLevenson et"all ([2009') argue that the observed 
small anisotropy in the mid-IR/X-ray correlation (see 
Sect. [S]) is qualitatively in agreement with a clumpy 
torus. Here we aim at being more quantitative and 
show consistency of the observed anisotropy with torus 
orientation. 

In Fig. [2]we overplot the observed ratio with predicted 
ratios of a clumpy torus model (red dashed line). As 
mean torus inclination in quasars we assumed 39° while 
the mean radio galaxy inclination is set to 75°. This 
corresponds to a mean opening angle of the torus of 60° 
or a type 1/type 2 ratio of 1:1, which is consistent with 
the number statistics of our complete and isotropic sam- 
ple (jLeipski et al..,2010,) . Fig. [2| shows that the model 



is following the continuum anisotropy curve quite well. 
There is, however, a slight deviation within the silicate 
feature where in the center of the feature the model 
curve predicts slightly lower anisotropy than shown by 
the observed curve. This may be indicative of addi- 
tional, off-torus obscuration in some objects (e.g. from 
the host galaxy), which effectively deepens silicate ab- 
sorption features in radio galaxies (the denominator of 
the plot in Fig. [2]) and, to a lesser degree, changes the 
spectral slope. If this happens in individual objects, the 
average curve and scatter will tend to be slightly more 
anisotropic. Some local examples of host obscuration are 
Centaurus A, NGC 5506, or the nucleus of the Circinus 
galaxy where host-galactic dust lanes are projected onto 
the nucleus producing deep silicate features. In fact the 
silicate absorption feature in the average radio galaxy 
SED seems to be deeper than in typic al local S eyfcrt 
2 galaxies without host obscuration ( Honig et al.|[2010,) . 
Note that this requires a statistical alignment of cold host 
dust with that of the torus (see Sect. I4.1|) . 

The model used for reproducing the ratio uses a dust 
cloud distribution with radial power law oc and 
5 clouds along the line-of-sight in equat orial direction 
(for details see iHonig fc Kishimotol I2010D . In compar- 
ison to models for SEDs and mid-IR interferometry of 
local Seyfert galaxies, these parameters suggest only a 
slightly more centr ally condensed and transparent torus 
(jHonig et al.ll2010() . while the half-opening angle may be 
wider (60° instead of 45°). In fact a range of torus 
model parameters satisfies the observed type 1/type 2 
anisotropy spectrum within the error bars of the sam- 
ple (e.g. various steeper and shallower dust distribu- 
tions). From Bayesian inference analysis we found that 
the torus model parameters are generally poorly con- 
strained (broad posterior distributions for individual pa- 
rameters). A weakness is certainly that modeling the 
flux ratio is not very constraining for model parameters 
since it does only take relative fluxes into account, while 
absolute fluxes (e.g. actual silicate strength of type Is or 
type 2s) are not included. On the other hand, what the 
modeling shows is that the observed small ratios at long 
wavelengths and the change of anisotropy from the near- 
to the mid-IR are in agreement with expectations from 
clumpy torus models without fine-tuning parameters. 

In summary, the torus model seems to reproduce the 
observed anisotropy reasonably well over most of the 
wavelength range, while single extinction laws result in 
much worse fits. The model parameters used in the 
clumpy torus model fit are reasonable in comparison to 
fits to local AGN. It demonstrates plausibility of the sce- 
nario that the anisotropy is related to torus orientation. 
This strongly suggests that the observed anisotropy is a 
measure for the intrinsic anisotropy of luminous type 1 
and type 2 AGN at z 1.2. 

5. OTHER ANISOTROPY ESTIMATES 

5.1. Comparison to the mid-IR/X-ray correlation 

In Sect, m we argued that the observed IR anisotropy is 
probably refiecting the "intrinsic anisotropy" as caused 
by the dust torus. The isotropic and complete selection 
of the sample helps to minimize any biasing effects on 
the anisotropy. On the other hand, optically- and X- 
ray-selected samples often suffer from missing some of 
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the most obscured objects. Moreover, since the ratio is 
> 1 at 15 /xm, flux-hmited mid-IR selected samples are 
potentially biased towards type 1 AGN as well. Thus, 
even using 12 /im-selected AGN samples will slightly suf- 
fer from the assumption of isotropy. 

One popular way of comparing type Is and type 2s is 
the correlation between X-ray and mid-IR luminosity. In 
case the X-luminosity is emitted isotropically (and traces 
the dust-heating emission), and the mid-IR emission is 
radiated very anisotropically, this correlation is expected 
to be different for type 1 and type 2 AGN. However , using 
high spatial resolution observations iGandhi et all ()2009n 
showed that local Seyfert 1 and Seyfert 2 galaxies, up to 
column densities of few 10^^ cm~^, essentially follow the 
same correlation. A conservative estimate suggests that 
at 12 /xm the difference between both types is smaller 
than a factor of 3. Despite including some "mildly" 
Compton-thick objects in this study, the most obscured 
objects are still missing due to the lack of intrinsic X-ray 
data. This essentially makes anisotropy estimates from 
the mid-IR/X-ray-correlation a lower limit on the "true" 
anisotropy. Nevertheless this result is fully consistent 
with our finding for powerful z ~ 1.2 AGN. At 12 /im 
we find a type 1/type 2 ratio of 1.65 ± 0.26. Assuming 
that X-ray selection misses the highest-inclination ob- 
jects our result would predict that the anisotropy in the 
mid-IR/X-ray correlation is < 1.6, at least for powerful 
AGN as presented here. 

5.2. IR anisotropy estimation is Seyfert galaxies 

iBuchanan et al.1 (|2006[ ) used a sample of local Seyfert 
galaxies and compared the average 5-35 fim SED of type 
Is and type 2s, scaled to their respective 8.4 GHz emis- 
sion. However since the sources have been selected ac- 
cording to a flux li mit at 12 /xm, the samp le cannot be 
considered isotropic. IBuchanan et al.l (|2006l ) report gen- 
erally higher fluxes for type Seyfert 1 AGN as compared 
to Seyfert 2s. The anisotropy decreases from a factor of 
about 8 to ~2.5 from 5 to 8 /im. This is significantly 
larger than what we find for our isotropically selected 
radio-loud sample. At longer wavel e ngths the Seyfert 
1 / Seyfert 2 ratio of IBuchanan et al.l ()2006l ) flattens to 
about a factor of 2-3 (with no convergence to unity as ex- 
pected at long wavelengths) which is, again, larger than 
what we found. 

The discrep a ncy between our results and 
IBuchanan et al.l (|2006f ) may be either due to (1) 
the different selection criteria chosen, (2) a difference be- 
tween low and high luminosity AGN, or (3) a difference 
between radio-quiet and radi o-loud AGN. While the 
data in IBuchanan et al.l ()2006f ) has not been corrected 
for host galaxy, the contribution by starformation should 
not be significant given the lack of PAH features in the 
average spectrum. If a host-correction were applied, it 
would predominantly affect type 2 AGN, thus making 
the anisotropy even larger. Furthermore, if the mid-IR 
selection had any effect, then the sample would miss 
out AGN at highest obscuration, so that the real type 
1/type 2 anisotropy would again be larger. In conclu- 
sio n, possible selection effec ts and host contamination in 
the IBuchanan et al.l ()2006l ) study would tend to result 
in an underestimated anisotropy, making the difference 
to our findings even stronger. 

It is well possible that the difference in anisotropy 



between our high-luminosity radio-loud sample and the 
low-luminosity radio-quiet sample is real. This would im- 
ply that either luminosity or radio power are the drivers 
for the observed characteristics. Radio jets are highly col- 
limated so that any influence of the jet can be expected 
perpendicular to the torus plane and just in a very small 
solid angle. It is, therefore, more reasonable to assume 
that the higher AGN luminosity would cause the lower 
anisotropy than the jet. The classical receding torus pic- 
ture changes thejjpcning angle of the torus for higher 
luminosity (e.g. Lawrence 1991; Simpson 2005). This 
mainly affects the relative number of type Is and type 2s 
in a sample. The relatively high fraction of about 50% 
unobscured AGN in our sample would support this sce- 
nario. However, to change the anisotropy between both 
types, it would be necessary to also change the obscura- 
tion properties (i.e. type 2s must on average look more 
like type Is). Such a scenario of "radiation-limited ob - 
scuration" has been proposed bv lHonig fc BeckertI (I2007D 
and supported by observations of lTreister et al.l ( )2009[ ). 
In this case about the same effect is expected for radio- 
quiet and radio-loud AGN. 

6. SUMMARY AND CONCLUSIONS 

We use the sample of qu asars and rad i o gala xies at 
z 1.2 recently presented in lLeipski et al.l (|2010[ ). Since 
the sample was selected isotropically, it should cover all 
torus inclination angles (weighted by solid angle). For 
these objects an infrared 1-17/im restframe SED has 
been constructed. Average SEDs were calculated for 
the quasar (= type 1 AGN) and radio galaxy (= type 
2 AGN) samples, respectively, in order to study the in- 
trinsic anisotropy of the IR emission of the dust torus. 

It is shown that the ratio between type 1 and type 
2 AGN in our parameter space of very luminous radio 
galaxies, the value gradually decreases from 20 to 2 at 
wavelengths 2 to 8 /im. Within the 10 /im silicate feature 
the ratio raises slightly. At longer wavelength the mid- 
IR emission becomes more isotropic. The intrinsic ratio 
between our type 1 and type 2 AGN is 1.44 ± 0.17 at 
15 /im. When using IR-selected flux-limited samples this 
anisotropy has to be taken into account. 

By analyzing the silicate feature in the sample av- 
erages we find the well-established "shift" of the cen- 
tral peak of the silicate emission feature with respect to 
the center of the absorption feature. The resulting cen- 
tral wavelengths at 9.6 ± 0.1 /im for the absorption and 
10.6 ± 0.1 /im for the emiss ion feature are in a.greement 
with previous reports (e.g. iSturm et al.ll2005l : iShi et al] 
I2OOI) . 

We discussed our results in the frame of other 
anisotropy estimators. Our findings are consistent with 
upper limits derived from the X-ray /mid-IR co rrelation 
of local Seyfert galaxies ([Gandhi et al.l 12009( 1. Some 
discrepancy exis t s wit h respect to a similar study of 
IBuchanan et al.l (|2006D for nearby Seyferts. If real it 
would imply that nearby, radio-quiet lower-luminosity 
AGN show a higher degree of anisotropy in the IR 
than higher luminosity, radio-loud sources. This may 
be explained by a receding torus model with luminosity- 
dependent obscuration. We also show that the overall 
relatively small degree of anisotropy is consistent with 
the torus being clumpy rather than smooth. Our clumpy 
torus model reproduces the observed type 1 / type 2 ratio 
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reasonably well. 
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